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K' < 21.0 sample is consistent with previous observations, and for the first time 
we derive N(z) down to K' = 24.0. After taking account of the dust extinction 
and selection effects of the sample, the observed N(z) distribution is well described 
with the PLE model, while the hierarchical galaxy formation model shows an 
apparent deficiency of galaxies especially at z > 2. The photometric redshift 
and the best-fit SED model evaluations allow the derivation of the rest-frame 
K', B, and UV(2000A)-band luminosity functions and their evolutions. The rest 
i^'-band luminosity function shows almost no evolution up to z = 3, while the 
rest B luminosity function shows mild luminosity evolution, and the rest UV 
luminosity function shows strong luminosity evolution. These findings seem to 
be qualitatively in favor of the PLE model. No evolution in the rest X'-band 
luminosity function can also be consistent with the hierarchical galaxy formation 
model if Mf Lk decreases with redshift, that is, if the number density of fT'-band 
luminosity- selected galaxies does not significantly change with redshift while the 
number density of stellar mass-selected galaxies decreases with look-back time. 
This trend corresponds to the evolution of the rest UV(2000A) — K' color, which 
gets bluer with increasing redshift. We also found that more massive galaxies are 
redder in this rest-frame color in every epoch. The rest-UV luminosity function 
of our i^'-selected galaxies shows a much shallower faint end slope at z = 3 than 
that of previous estimations for rest-UV selected Lyman break galaxies. As a 
consequence, the contribution to the global star formation rate of our i^'-selected 
galaxies is about 42% of that derived from the integration of luminosity function 
of Lyman break galaxies at z — 3. This result suggests that a large fraction of 
the star formation rate density at z > 1.5 comes from the contribution from the 
faint (M 2 ooo > —20) blue galaxy population at high redshift universe that have 
not yet obviously been identified. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
luminosity function, mass function 

1. Introduction 

The large galaxy sample with redshift determined to faint magnitudes allows us to 
understand the evolution of galaxies. Previous deep spectroscopic surveys have shown several 
evolutionary trends for galaxies at < z < 1, especially for the evolution of luminosity 
functions (LFs), which is one of the fundamental observational tools to trace the galaxy 
evolution. Lilly et al. (1995) have shown that the LF for red galaxies does not evolve 
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significantly over < z < 1, which infers that red massive galaxies must have already been 
assembled by z = 1. On the other hand, their blue subsample is evolving, which indicates 
that they are formed later than the red subsample. Cowie et al. (1996) have shown the 
evidence of the "downsizing" trend in which more massive galaxies appear to be forming at 
higher redshifts. 

Although it is important to see these evolutionary trends in the earlier epoch beyond 
z — 1, it is difficult to construct a large galaxy sample in the familiar optical wavelength 
except for the distinct high z star-forming populations of Lyman break galaxies and/or Lya 
emitters at z > 3 (Steidel et al. 1999; Ouchi et al. 2003, Paper II). There are some possible 
interpretations of such luminous high z galaxies such as a core of today's massive L > L* 
galaxies (Steidel et al. 1996) or as a low-mass starbursting building blocks of the more massive 
galaxies seen today (Lowenthal et al. 1997). However, it remains uncertain how these z > 3 
galaxies are related to present-day galaxies. 

One of the explicit difficulties in tracing systematically the evolution of galaxies is that 
a single observed frame does not correspond to a single rest frame over wide redshift range, 
say < z < 3. The rest-optical frame in which we have seen the evolution trend at < z < 1 
shifts to the NIR region at z > 1.5. An observation in wide wavelength range is undoubtedly 
required to know the SEDs of galaxies and follow their evolution based on the same rest 
frame in which we can draw almost similar stellar contents at any redshifts. 

There are two approaches to probe the rest-frame luminosity of galaxies. One is to 
change the observed band accordingly and follow the redshifting rest frame. This straight- 
forward approach was taken by many studies. Cowie, Songaila, & Barger (1999) traced the 
rest-S band up to z = 1 by UBVI multiband observation. Furusawa et al. (2003, Paper 
IV) used multicolor BVRi'z' bands and traced the rest B band LF up to z = 1.25. NIR 
photometry for Lyman break galaxies enables them to derive their rest V^-band luminosity 
at z = 3 (Shapley et al. 2001). This approach gives a reliable way to trace the rest frame 
directly, though the study is restricted to a single rest band. 

The other approach is to utilize the best-fit SED templates determined on photometric 
redshift method (Sawicki, Lin & Yee 1997; Rudnick et al. 2001). As have been already 
mentioned by many authors, the photometric redshift technique has less precision for redshift 
estimates than the spectroscopic redshift. However, spectroscopic redshifts are practical only 
for bright objects. Photometric redshifts can be extensively applied to fainter objects with 
a fair degree of confidence. They are also of great use in constructing a statisticaly large 
sample in the redshift range between z = 1.5 and z = 2.5, where redshifts can be identified 
only by NIR spectroscopy. It works predominantly as a statistical technique rather than as 
a constraint on parameters that determine the SED of an individual object. Studies based 
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on the photometric redshifts and those based on the spectroscopic redshifts complement 
each other. If the SED of each sample galaxy in wide wavelength is given, the rest-frame 
luminosity of the galaxy can be estimated. This method has the advantage to deriving the 
rest-frame luminosity at any redshift and at any wavelength, though the estimated luminosity 
inherits an uncertainty in determining the photometric redshift. 

In this study, we have constructed a very deep K' band-selected sample of galaxies. 
We then adopt the second approach described above with photometric redshifts derived 
from multicolor BVRIz'JK' deep images and investigate the evolution of the LFs in the 
three (K', B, and UV) rest-frame bands out to z — 3.5 where the i^'-band still traces the 
rest wavelength longer than an age-sensitive spectral break at 4000A. NIR band-selected 
samples may provide significant advantages over optical band-selected samples in studying 
galaxy evolution due to smaller extinction by dust and less type-dependent k— corrections 
(Mannucci et al. 2001) at these wavelengths. Moreover, NIR selection provides samples that 
are not biased towards star-forming galaxies and allows estimation of the mass of galaxies 
over a wide range of redshift (Saracco et al. 2001). NIR photometry in general improves the 
estimate of photometric redshift (e.g., Bolzonella, Miralles, & Pello 2000). 

The Subaru Deep Field (SDF) was imaged quite deeply at near-infrared with wide field 
near infrared camera CISCO and yielded some of the faintest galaxies ever observed, down 
to magnitudes of J = 25.5 and K 1 = 24.5 (Maihara et al. 2001, Paper I), which is much 
deeper than most of previous i^'-band surveys (e.g., Cowie et al. 1996; Barger et al. 1999; 
Fontana et al. 2000; Rudnick et al. 2001). The field was well selected, taking care to avoid 
large Galactic extinction and nearby clusters of galaxies and to get smaller airmass than the 
HDF at Mauna Kea (Paper I). The following optical imaging enabled us to estimate the 
photometric redshifts of these faint objects and offered the discovery of faint high z object 
candidates. In the SDF, we have already evaluated the number counts in J and K' bands, the 
near-infrared color distribution, and contribution to extragalactic background light (Paper I, 
see also Totani et al. 2001a). Totani et al. (2001b) have investigated a detailed comparison of 
the number counts, colors, and size distribution for the SDF galaxies with theoretical galaxy 
formation models. Nagashima et al. (2003) have also analyzed the SDF galaxies by using 
a semianalytic model of galaxy formation based on the CDM model. In this paper, we will 
compare the redshift distribution with their models. 

This paper will be organized as follows. We present the process to construct the K'- 
limited sample of SDF galaxies and estimating photometric redshifts as well as their spec- 
troscopic calibration in §2. In §3, we show the redshift distribution of our sample and 
compare it with galaxy formation models. We discuss the evolution of LFs on rest K', B, 
and UV(2000A) bands in §4 and the evolution of rest UV— K' color in §5. We will also 
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describe in §6 the evaluation of star formation rate density (SFRD) derived from our LF 
estimate. We give a summary in §7. Throughout this paper we assume a flat, matter domi- 
nated universe (Q m = 1) with H = 75kms _1 Mpc~ 1 to compare with other major previous 
observations. 



2. K' band-selected Galaxy Sample 
2.1. Data 

We have used a deep multicolor data set of BVRIz' JK' images taken with the Subaru 
telescope on the Subaru Deep Field (SDF). The SDF is centered on a = 13 /l 24 m 21. s 38; 5 = 
+27°29'23."0(J2000). The J and K' imagings and data reductions are described in the 
previous paper (Paper I). The optical follow-up imaging observations were carried out during 
the commissioning runs of the FOCAS (Kashikawa et al. 2002) and Suprime-Cam (Miyazaki 
et al. 1998) in 2001 February to June. In this study, we used the Johnson B-, Kron-Cousins 
R- and /-band images taken with FOCAS and Johnson V- and SDSS z'- band images taken 
with Suprime-Cam. The observation and data reduction of the V and z' data, which are one 
of the wide field imaging data for the SDF, will be described in coming paper (Furusawa et 
al. 2003, Paper IV). Here we briefly describe the imaging data of B, R, and I taken with 
FOCAS. 

FOCAS has a field of view as large as 6'0, which could entirely cover not only the SDF 
2' x 2' itself but also its flanking fields. The SDF has flanking fields 25 arcmin 2 in total size 
surrounding the SDF, though their limiting magnitude in K' (K' < 20.5) is fairly shallower 
than that of the central field. FOCAS has a couple of 2048 x 4096 pixel CCD chips with 
the pixel scale 0.103 arcsec/pix. The typical seeing size was 0.6 — 0.8 arcsec. The total 
integration time was 2.5hr, 4hr, and lhr for B, R, and / bands respectively. 

The B, R, and / images were reduced in a standard manner. After flat-fielding, we 
have corrected the image distortion caused by the optics of FOCAS and the telescope. The 
undistorted images have accuracy of 1.0 pix rms. This distortion pattern has been confirmed 
to be independent of the flexure of the telescope and the instrument. The image shift caused 
by the flexure was about 0.6 pix within the maximum exposure time of 30min. Relative 
sensitivity of two CCDs was corrected based on the difference of sky brightness assuming 
that the sky is uniform over the field of view. Actually, the relative ratios of sky brightness of 
two CCDs on the same exposure were stable with 1% fluctuation. The sky was determined by 
fitting the unvignetted area of each frame with two-dimensional 3rd or 4th order polynomials 
and subtracted. The relative shifts and rotations among dithered frames were measured and 
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fixed using the positions of several point sources. Then they were stacked each other using 
the IMCOMBINE task in IRAF. Magnitude zero points were calibrated to the Johnson- 
Cousins system using Landolt photometric standard stars. The accuracy of the zero points 
is estimated to be in the range 0.03 — 0.08mag. The magnitudes have also been corrected 
for the Galactic absorption, using E(B — V) — 0.018 derived from Schlegel, Flinkbeiner & 
Davis (1998), yielding A B = 0.074mag, A v = 0.056 mag, A R = 0.046 mag, Ai = 0.033 mag, 
A z , = 0.025mag, Aj = 0.015 mag, and A K , = 0.006 mag. 

The summary of optical imaging data used in this study is presented in Table 1. The 
5<T limiting magnitudes in Vega system within 2."0 aperture for these images are B < 28.08, 
V < 26.28, R < 27.88, / < 25.93, z' < 24.27, J < 25.11, and K' < 23.73. 

Among seven images, the J-band image has the smallest area, 2' x 2', on which the 
deep i^'-band imaging has carried out. Therefore we first rescaled all the data to the same 
pixel scale (0."103 /pix) as the B, R, and / data, corrected relative rotation based on the 
positions of several bright point sources, and extracted the 1/90 x 1/97 region common to 
the J- and K'-band images for further analysis. The error of this image translation was 
below 0.02 arcsec. 



2.2. Object Detection and Photometry 

We extracted objects using SExtractor (Bertin & Arnouts 1996). The objects were 
extracted from the K'-band image as a reference image and photometric measurements 
were performed on the other image for the exact positions where sources were detected 
on the reference image. We choose a 4.0pix FWHM Gaussian convolution kernel in object 
detection on the reference if'-band image. We used a detection threshold of 24. 1 mag/ arcsec 2 
corresponding to the 1.5a above the background noise level in K' band, the same as in Paper 
I. 

To measure colors of detected objects, aperture photometry was carried out for seeing- 
corrected images, which were created by convolving Gaussian filters to match the seeing 
size to that of the V band, the largest among the seven images. At first we used a fixed 
aperture. However, it was difficult to choose an adequate aperture size because photometry 
with a large aperture is affected by close neighboring objects while a small aperture does not 
cover the same physical sizes for both nearby objects and distant objects at the same time. 
Finally we chose an aperture size adjusted for each object size but kept identical for all seven 
images. To choose the adequate aperture size in deriving colors of objects, we have tried 
many times with different aperture sizes and confirmed the consistency between the resulting 
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photometric redshifts and spectroscopic redshifts for calibration samples (see Section 2.3). 
We found that 2xFWHM (FWHM was measured on i^'-band image) aperture for each 
object is the optimum aperture size to derive correct photometric redshifts. In this way, we 
have measured total magnitudes ("MAG_BEST" magnitutudes from SExtractor), 2 x FWHM 
aperture magnitudes for all the seven images. It should be noted that the 2 x FWHM aperture 
magnitude is only used for SED fitting in deriving the photometric redshift, and further 
analysis is based on the total magnitude. The nine objects that have stellarity> 0.9 with 
K' < 23.0 were eliminated as Galactic stars from our galaxy sample. We have also removed 
12 objects whose photometry is strongly affected by nearby bright stars in one of the seven 
band images. 

Our final i^'-selected catalog contains 439 objects with the magnitude limit K' tot < 24.0. 
According to Paper I, the completeness is over 70% at this limiting magnitude. 

2.3. Photometric Redshifts 

We used the publicly available Hyperz code by Bolzonella, Miralles, & Pello (2000) to 
determine photometric redshifts for our seven-band BVRIz' JK' photometry. The Hyperz 
code derives a photometric redshift for each input multi-photometry set according to the 
standard SED fitting minimization procedure. The detail of the method as well as discussions 
about the accuracy of evaluated redshifts are presented in Bolzonella, Miralles, & Pello 
(2000), therefore we describe here briefly about our choices of major input parameters. We 
used five (Burst, E, Sa, Sc, and Im) synthetic galaxy SED templates of Bruzal & Chariot 
model (GISSEL98, Bruzal & Chariot 1993) to compare with the observed SED. The default 
SED templates applied in Hyperz assume the Miller & Scalo IMF from 0.1M Q to 125M , 
solar metallicity Z = Z & = 0.02, ages ranging from 10 6 yr to 2 x 10 10 yr, and Madau (1995)'s 
formula for continuum depression due to the intergalactic absorption blueward of Lyman 
a emission. The adopted SED templates are assumed for an instantaneous single-burst 
SFR model for "Burst" SED, and exponentially decaying SFRs with e-folding timescales of 
1, 3, 15 Gyr and infinity (i.e., constant SFR) to match the colors of present-day galaxies 
of different spectral types from "E" to "Im". Though three more additional spectral type 
templates can be applied in Hyperz, the accuracy for photometric redshift z p hot estimate was 
found to be unchanged when applying all the eight templates. We have chosen a Calzetti 
etal. (2000) 's extinction law with Ay ranging between 0.0 and 1.8. The step in redshift in 
searching the solutions between z = and z = 7 is Az = 0.1, though the resolution in the 
final minimization is 10 times better. The results are not significantly affected by the choice 
of smaller Az. We have also constrained on the age of the template to be less than the age of 
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the universe at the redshift under consideration. The transmission curve for a given filter set 
was defined based on the real measured curve for each instrument. An uncertainty of derived 
photometric redshift for each object was estimated from the x 2 probability distribution at 
the confidence level 90% derived in SED template matching, which will be shown later. As 
expected, the error becomes larger as K' magnitude goes fainter. 

We examined the accuracy of our determination of photometric redshifts, z p hot, against 
spectroscopic redshifts, z spec . We have measured spectroscopic redshifts for 76 galaxies in 5' x 
5' field centered on SDF. All these spectroscopic data were taken during the commissioning 
runs of FOCAS MOS mode in April and May, 2001. The spectra cover 4700 - 9000A with 
a pixel resolution of 2.8A under the 0."8 — l/'O seeing size and were reduced in a standard 
manner. 

In this 5' x 5' SDF flanking fields, J-band data are not available except for the central 2' x 
2' region and the if'-band image has shallower limiting magnitudes, K' < 20.5. Five galaxies 
of 76 galaxies belong to the SDF central field in which J-band photometry is available. We 
have confirmed an excellent agreement of the estimated z p hot between the photometric data 
sets with and without J-band for these five galaxies. The comparison of z p h ot with z spec for 
our spectroscopic sample is shown in Figure 1, in which we found good agreement, though 
we don't have a large enough spectroscopic galaxy sample at 1.5 < z. The mean error is 
Az = (z spec — z p hot) = —0.038 with a = 0.207. The large inaccuracy seen in z < 0.4 would be 
caused by the lack of U— band photometry which can detect a Balmer break at this redshift 
range. If objects with z < 0.4 are removed, the accuracy improves to a = 0.178. Note that 
we don't have any apparent catastrophic errors, though the spectroscopic sample itself is 
small. If five objects with Az > 0.4 are removed, the accuracy improves to a = 0.160. 



3. Redshift Distribution 

We describe in this section the derived redshift distributions of our sample. Figure 2 
shows the K'-z relation. The error bars of the determined z p h ot values are estimated from 
the x 2 probability distribution at confidence level 90% derived in SED template matching. 
Figure 3 shows the redshift distribution of our if' -selected sample. The distribution has a 
peak around z = 0.8 which is slightly larger than the N(z) peak around z = 0.5 — 0.6 derived 
by Fontana et al. (1999) or Rudnick et al. (2001) for K < 21 samples which include HDF-N,S 
fields. However, if we restrict our sample to brighter galaxies having K' < 21.0, the peak 
moves slightly toward a lower z at z = 0.6. The number density of galaxies at 1 < z < 1.5 
is 4.8 galaxies/arcmin 2 at K' < 21.0, which is consistent with the result of Fontana et al. 
(1999) of 4.4 galaxies/arcmin 2 at K' < 21.0. Thus, our N(z) is consistent with those of 
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Fontana et al. (1999) and Rudnick et al. (2001). 

We can use this overall redshift distribution of K' band-selected galaxies to compare 
the predictions of theoretical models of galaxy formation. The redshift distribution for the 
if-selected galaxy sample was first pointed out by Kauffmann & Chariot (1998, KC98) as an 
effective tool to discriminate amoung galaxy formation models. The hierarchical clustering 
model, in which massive galaxies were assembled much more recently through merging of 
small galaxies, predicts less galaxy abundance at high z (z > 1) than that expected from 
the pure luminosity evolution (PLE) model. The recent results of Fontana et al. (1999) and 
Rudnick et al. (2001) for the K' < 21 sample were both consistent with the hierarchical 
clustering model. Since our if' -selected sample is much deeper than Fontana et al. (1999) 
and Rudnick et al. (2001) by 3mag, the derived redshift distribution reflects both a high z 
(z > 2) population and intrinsically faint population in the local universe and is crucial and 
useful in reconciling observation with theoretical models. Actually, the fraction of galaxies 
at z > 1.6 detected in K' < 21 sample of Fontana et al. (1999) was only 9%. 

Following the same manner as KC98 and previous observations, Figure 4 demonstrates 
the normalized cumulative redshift histogram for our i^'-selected sample, as well as the the- 
oretical predictions of both the PLE and hierarchical clustering models. The thick solid 
histograms show the results of observation, while other lines show the model predictions. 
We have carried out Monte Carlo realizations to see the distortion that the error of our 
derived photometric redshift could cause on the resulted N(z). At first, we have created a 
"pseudo-galaxy-catalog" in which redshift of each galaxy was assigned to the photometric 
redshift solution plus a random error perturbed in the 90% range of confidence interval in the 
SED x 2 fitting procedure (same as the error bars denoted in Figure 2). In this process, mag- 
nitude of each galaxy was also assigned with a random error perturbed within the measured 
photometric error by SExtractor. We assumed here that the error in photometric redshifts is 
independent of the photometric error, though a large photometric error would cause a large 
error of photometric redshift in most cases. The redshift distribution is recalculated for this 
sample. We repeated the process of generating such pseudo-galaxy-catalogs 200 times and 
investigated the scatter of derived redshift distribution. The shaded regions in Figure 4 show 
the ±3<7 deviated counts estimated by this method. 

The thick dashed lines show the PLE model by Totani et al. (2001b). The model 
normalizes its number density of galaxies at z = by the B band luminosity function, and 
it probes the evolution backward in time. Galaxies are classified into six morphological 
types of E/S0, Sab, Sbc, Scd, Sdm, and the distinct population of dwarf elliptical galaxies. 
Their luminosity and SED evolutions are described by a standard galaxy evolution model of 
Arimoto & Yoshii (1987) and Arimoto, Yoshii & Takahara (1992). All galaxies are simply 
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assumed to be formed at a single redshift, zf = 5. A reasonable evaluation of dust extinction 
by a screen model, and selection effects such as cosmological dimming of surface brightness 
are taken into account. The details of the model are described in Yoshii (1993) and Totani 
& Yoshii (2000). Their selection criteria for object detection are completely the same as ours 
and those in Paper I. This PLE model, which takes account the two distict populations of 
giant elliptical galaxies and dwarf elliptical galaxies, was found to match extremely well with 
the observed K-band galaxy counts on the SDF beyond K > 22.5 where the PLE model 
with a single elliptical population could not explain the observed galaxy abundance. In their 
original model of Totani et al. (2001b) they have predicted the N(z) for the SDF based on an 
isophotal magnitude in which selection effects are easy to take into account. We recalculated 
the N(z) predicted by the model based on the total magnitude in this study. 

On the other hand, the thick dotted lines represent the semianalytical hierarchical galaxy 
formation model by Nagashima et al. (2003). This model is based on the structure formation 
theory in the CDM model, which corresponds to the hierarchical clustering of dark halos, 
and contains some important physical processes concerned with the evolution of the baryonic 
components such as gas cooling, star formation, supernova feedback and galaxy mergers. In 
this model, merging histories of dark halos are realized by using a Monte Carlo method 
assuming a power spectrum of initial density fluctuation in a Gaussian random field, and 
then, the above baryonic processes are followed in each dark halo. Thus star formation 
histories, luminosities and colors of individual galaxies are predicted with this model. The 
details of the model itself are presented in Nagashima et al. (2001) in which the model well 
reproduces the local luminosity function, the H I gas mass fraction, galaxy sizes and so on. 
The Nagashima et al. (2003) model employed in this study improves on the Nagashima et 
al. (2001) model by taking account almost the same selection effects as those in Totani et al. 
(2001b) and dust extinction, and can predict the number counts of galaxies consistent with 
that of observed on the SDF. Note that the model of Totani et al. (2001b) was calculated 
for a cosmology with Vt m = 0.2 and Q\ = 0.8 and that of Nagashima et al. (2003) was 
calculated with Vl m = 0.3 and VL\ = 0.7 as were those of their original models, whereas 
an overall distribution of observed photometric redshift does not depend on cosmology. In 
the brightest 19 < K' < 21 bin, we also overlay the model predictions of the PLE and 
hierarchical clustering models of KC98 for comparison. 

It should be noted that the hierarchical clustering model apparently shows a strong 
deficiency compared with the observed N(z), especially at high redshifts (z > 2). In the 22 < 
K' < 23 bin, which contains 119 sample galaxies and still attains the 100% completeness, 
the PLE model agrees with the observed count better than hierarchical galaxy formation 
model does. The PLE model can account for the abundance of the high-z (z > 2) galaxies. 
On the other hand, its underprediction of the high-z galaxy abundance as compared with 
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that of the PLE model of KC98 is attributed to the dust evaluation in the model. 

We performed a two-sample form of the Kolmogorov-Smirnov test (K-S test) to see the 
statistical agreement between observation and models. The chance probability, as well as the 
X 2 values (in parentheses) are listed in Table 2. The K-S test accepts at the 1% significance 
level the hypothesis that the observation and the PLE model in the 19 < K' < 23 bins 
are drawn from the same distribution, while the hierarchical galaxy formation model agrees 
with observation only in the bright 19 < K' < 22 bins. In the faintest 23 < K' < 24 bin 
both the model predictions are found to be different from the observation. In this bin, the 
observed count shows a high z tail which may be an effect of having more galaxies at high 
z and fewer galaxies at low z than in the models. Though the error-shaded regions seem 
to be distributed below the observed counts, it is partly because the pseudo-galaxy catalog 
lacks the objects whose reassigned magnitudes are beyond the limiting magnitude K' = 24.0. 
Note that the completeness of the sample is 100% at K' = 23.5, though 70% at K' = 24.0 
(Paper I). In the galaxy count, the model prediction begins to show a slight deviation from 
the observation at this faint magnitude (Totani et al. 2001b; Nagashima et al. 2003). The 
discrepancy in this faintest bin is presumably caused by the large photometric errors which 
give incorrect solutions of photometric redshifts. More detailed refinements of the model 
may also be required. 

As noted above, previous works on the KC98 test such as Fontana et al. (1999) and 
Rudnick et al. (2001) found that their N(z) is well described with the hierarchical clustering 
model prediction and that the PLE model overpredicts the abundance at all redshifts. Our 
observed N(z) itself is consistent with those observed by Fontana et al. (1999) and Rudnick 
et al. (2001) at least at K' < 21, and the KC98 test with our N(z) would give the same 
conclusion if we adopt the model predictions of KC98. In this study, we introduce more 
realistic models of PLE and hierarchical galaxy formation that take account of both the se- 
lection effects of the sample and dust extinction, and found that the PLE model is favored by 
the observed N(z), while the hierarchical galaxy formation model prediction is substantially 
below the observed N(z) at high z. Cimatti et al. (2002) concluded that their spectroscopic 
redshift distribution at K s < 20 of the K20 survey is consistent with the PLE models that 
take into account the selection effects and dust extinction. It might suggest that we need 
some additional physical processes in the CDM model. We note, however, that the result is 
obtained on the very small field of the sky, therefore it should be confirmed by a future deep 
and wide-field ^-selected sample. 
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4. Evolution of Rest-frame Luminosity Functions 

Using the estimated photometric redshift catalog, we now investigate the galaxy pop- 
ulation and its evolution by means of luminosity functions. Our wide wavelength coverage 
and SED estimate for each sample galaxy allow us to derive LFs in the three rest-frame K', 
B, and UV(2000A) bands. We will investigate the evolution of rest-frame LFs at z < 3.5, 
where K'-hand traces the rest-frame wavelength longer than 4000A. Though it is quite in- 
teresting to see the difference of rest-frame LFs for a blue or red subsample (e.g., Lilly et al. 
1995) or a spectral type (e.g., Folkes et al. 1999) subsample, our sample is not so large as 
to be subdivided into small subsamples with enough statistical confidence. Such a detailed 
analysis would be possible in the future when a large NIR imaging survey is carried out. 

4.1. Rest K' Frame Luminosity Function 

First, we derive the rest K'-band luminosity functions. For each galaxy we have deter- 
mined a rest-frame absolute K' magnitude, Mk> as 

M K , = m K > - 51og(d L /10pc) - k(z), (1) 

where d^ is the luminosity distance and uik 1 is the apparent K' magnitude. The 
A;— correction k(z) can be derived from the best-fit SED template for each galaxy. 

To determine the LFs we use a traditional l/V max method (Felten 1977), 

where dV/dz is the comoving differential volume. The limits z m i n and z max are defined 
as the smallest and the largest redshifts, respectively, that a given galaxy could have and 
still be included in a sample (see detail in Ellis et al. 1996). The sum is carried out over 
all galaxies in the sample lying within the specified range of redshift and magnitude. The 
volume is based on an effective solid angle, da;, of the survey of 3.739 arcmin 2 . 

The rest K' frame LFs for four redshift intervals at 0.6 < z < 3.5 are shown in Figure 5. 
We have neither enough sample galaxies to derive a reliable LF nor good accuracy in redshift 
identification at z < 0.6, so we omit this redshift interval. Instead, we refer the local K-hond 
LF derived by Loveday (2000) for the Stromlo-APM Redshift Survey, shown as a dotted line 
on each frame in Figure 5. Although Loveday's sample is smaller than that used in another 
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local K-band LF of 2dF survey (Cole et al. 2001), it extends to a magnitude as faint as 
M K = —16, with smaller statistical error, which is adequate to compare with our sample. 
The error bars show the statistical uncertainties evaluated from a bootstrap resampling 
method with a series of 500 Monte Carlo simulations. We fit the LF by a usual Schechter 
parameterization and list the best-fit parameters in Table 3. 

The derived K' LFs at 0.6 < z < 1.0 differ slightly from the local one, having hints of a 
steeper faint-end slope and brighter bright end, though the overall shape of the LF seem to 
agree with that of local i^'-band LF taking account the normalization difference. The most 
striking result perceived from this figure is that the i^'-band LF remains almost unchanged 
from z = 0.6(0.0) to z = 2.5 and shows decline at the bright end at 2.5 < z < 3.5. We have 
evaluated the robustness of our resulted LFs against the errors of our photometry and the 
photometric redshift estimate by the Monte Carlo realizations, which is as the same as those 
we performed in §3. The evaluated la scatter of derived LF parameters, i.e., a, M*, and <p* 
are denoted in Table 3. We found that the effects of photometric and redshift errors are small 
with negligible effect on a, and at most 0.5mag in M*. The i^'-band LF at 2.5 < z < 3.5 
also has signs of a steeper faint end slope than that at lower z ones, which is expected in 
the context of the hierarchical clustering model. However, the trend may not be statistically 
significant because of the lack of enough sample objects in faint [Mk 1 > — 21) bins. 

Cowie et al. (1996) have derived the rest i^'-band LF for 0.6 < z < 1.0 galaxies. It 
remains almost unchanged from that for < z < 0.2. Drory et al. (2001) also showed no 
significant evolution of rest JT'-band LF for z < 1.2. Lilly et al. (1995) have also noted that 
rest-frame red luminosity galaxies have a weak evolution up to z = 1. This trend has been 
found to continue up to z < 3 in this study. The almost no evolution in the rest K'-band LF 
implies that the old stellar population in galaxies was already in place by z = 2 — 3 and that 
their stellar mass did not change significantly at < z < 3. It seems that the present result 
is inconsistent with the hierarchical clustering model. Papovich, Dickinson, & Ferguson 
(2002) estimated the stellar masses for the HDF-N galaxies by using population-synthesis 
models and found that massive present-day galaxies were not fully assembled by z ~ 2.5 
as predicted by the hierarchical galaxy formation model, which might be also inconsistent 
with ours. However, Drory et al. (2001) have shown that the M/L K slightly decreases with 
redshift at z < 1.2. If this trend continues up to z ~ 2.5, all the findings are consistent 
with the fact that the number density of K' band luminosity-selected galaxies does not 
remarkably change with redshift, while the number density of stellar mass-selected galaxies 
decreases with look-back time. 

It should be noted that we are assuming here a rather large extrapolation along the SED 
template from the observed i^'-band to the rest if'-band especially at 2.5 < z < 3.5. The 
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reliability of this extrapolation depends on some appropriate assumptions in constructing 
the SED templates, based on our fundamental knowledge of star formation histories in a 
galaxy. Recently, Shapley et al. (2001) and Papovich, Dickinson, & Ferguson (2002) have 
attempted to derive the stellar masses for individual sample of Lyman break galaxies at 
z ~ 3. In these studies, their UV to near-IR photometries were fitted with a series of stellar 
population synthesis models, which is basically the same SED-fitting method as photometric 
redshift except that they considered carefully almost full allowable range of various input 
parameters in constructing SED templates. This constraint method can be applied only 
to the samples whose redshifts are reliably determined by spectroscopy. Although some 
interesting trends have been found among parameters of star formation histories, they found 
only loose constraints on each of the parameters including stellar masses. This means that 
the stellar mass measurement even at observed K' band at z ~ 3 has some uncertainty. 
The most significant uncertainty factor in deriving stellar mass is a second stellar population 
component (Papovich, Dickinson, & Ferguson 2002) of old stars that have formed in a 
burst that occurred before ongoing active star formation. The stellar masses with this two- 
component model are estimated to be 3 — 8 times of that of the single-component models. 
While this old population, if any, actually influences on the estimate of stellar masses, its 
allowable flux contribution is supposed to be small because of its high mass-to-light ratio. 
For example in Figure 19 of Papovich, Dickinson, & Ferguson (2002), which shows the effect 
of a second old stellar population on the galaxy SED, the second component increases the 
inferred rest if' -band luminosity by a factor of less than ~ 1.4, which corresponds to ~ 0.37 
mag brightening. As a consequence, our extrapolation of SED and derived rest i^'-band 
magnitudes at high z in this study could have uncertainty to the same extent. 

In the following subsections, we will see the evolution of rest-frame LFs of bluer rest 
bands that trace young stellar population and are sensitive to star formation activity. 

4.2. Rest B Frame Luminosity Function 

We have also estimated the rest .B-band luminosity functions with almost the same 
kind of K' LF. Six objects out of 439 in the i^'-selected sample were not detected in the 
.B-band image. We have derived their inferred -B-band magnitudes from their best-fit SED 
templates, convolved with a .B-band response curve. The consistency between photometric 
-B-band magnitudes and these inferred B magnitudes from their best-fit SEDs were checked 
with sample galaxies except these six 5-dropout galaxies, and their agreement was found to 
be within 0.56 mag rms. The inferred B magnitudes of the six B dropout galaxies were all 
beyond the la limiting magnitude 29.76 of the -B-band image. 
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The rest £?-band LFs at 0.6 < z < 3.5 are shown in Figure 6, and the best-fit Schechter 
parameters are listed in Table 3. Superposed dotted lines are the local 5-band LFs from 
the SDSS (Blanton et al. 2001) converted from g*-band. In the 0.6 < z < 1.0 bin panel, we 
also superposed as a dashed line the 0.75 < z < 1.0 rest 5-band LF derived in the CFRS 
(Lilly et al. 1995) for the Q m — 1, Q\ — 0, and h = 0.75 case. Although CFRS samples only 
the bright end of the LF, the agreement between CFRS and our result for 0.6 < z < 1.0 
is quite good, both show the brightening by 0.5 mag from the local rest B LF. Note that 
CFRS is the /-band magnitude-limited sample and ours is if '-band limited, though both of 
these bands trace above 4000 A at z < 1.0. 

In the 2.5 < z < 3.5 panel of Figure 6, the dashed line represents the rest V-band LF 
of Lyman-break galaxies at z ~ 3 by Shapley et al. (2001). We have converted their rest 
V-band magnitude to our rest 5-band magnitude by using the rough estimate of average 
color of observed (H — K) ~ 0.40, which approximately corresponds to the rest (B — V), 
at z ~ 3 derived on HDF-N galaxies by Papovich, Dickinson, & Ferguson (2002). We have 
also reestimated the M v and <fi* values assuming a cosmology with Q m = 1.0 and h = 0.7. 
Though band conversion has some uncertainty (~ ±0.5), our rest B LF at z ~ 3 is in good 
agreement with that of Shapley et al. (2001) which is part of the rest UV spectroscopically 
selected sample of Steidel et al. (1999). The agreement of LFs at z ~ 3 of our if' -selected and 
their optically selected samples means that there are very few, if any, red candidates of old 
non-star-forming galaxies that would not be detected in optical bands but in near-infrared 
bands, which is also suggested by Shapley et al. (2001) and Papovich, Dickinson, & Ferguson 
(2002). On the contrary, we can conclude that our if '-selected sample does not have a large 
deficit of young blue star-forming galaxies at least in brighter magnitudes. 

In contrast to the rest K' LF, the rest B LF shows brightning with epoch up to z < 2.5. 
The bright end of the LF appears to continue brightening by 0.4mag from 0.6 < z < 1.0 
to 1.0 < z < 1.5, and by a further 0.4mag from 1.0 < z < 1.5 to 1.5 < z < 2.5, keeping 
the faint-end slope constant, though its brightening degree is within the M* error caused by 
photometric and redshift errors. The bright-end amplitude of the rest B LF at 1.5 < z < 2.5 
clearly exceeds over the local 5-band LF, while it declines at 2.5 < z < 3.5 almost to the 
local amplitude. 

This brightening evolution in rest B LF is also seen in Sawicki, Lin & Yee (1997) for 
the HDF F814Wab < 27 sample. They have found not only the brightening but also the 
steepening of the faint end slope with increasing redshift to z — 3. The difference in faint- 
end slope is presumably caused by the difference of the observed band used to construct the 
magnitude-limited sample, which will be discussed later. 

In the rest B band, we have found a strong luminosity evolution at z < 3 which was not 
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seen in the rest K' band. In next subsection, we further investigate this trend in the bluer 
UV band, which directly relates to the star formation rate of a galaxy. 

4.3. Rest UV 2000A Frame Luminosity Function 

We derive the rest UV LF at 2000A for which Sullivan et al. (2000) have derived the 
local LF for an ultraviolet selected sample. We have evaluated photometric redshifts, as 
well as best-fit SEDs, for all the galaxies in our K'-selected sample. Therefore we can 
derive absolute magnitudes of any rest wavelength in the AB system, although the derived 
magnitudes would suffer directly the error of the photometric redshift evaluation. Here we 
derive the absolute magnitude, M 2 ooo as follows, 

M 200 o = m AB + 2.51og(l + z) - 51og(d L (z)/10pc), (3) 

where m^ is the observed redshifted magnitude at 2000(1 + z)A and d^z) is the 
luminosity distance. It is to be noted that the derived rest UV luminosity is uncorrected 
for dust extinction, which is estimated to be in the range of Ai 70 o = 0.0 — 4.0mag at z ~ 3 
(Papovich, Dickinson, & Ferguson 2002). 

The rest 2000A LFs at 0.6 < z < 3.5 are represented in Figure 7 and best-fit Schechter 
parameters are again listed in Table 3. Superposed dotted lines are the local rest 2000A LF 
(Sullivan et al. 2000) for h = 0.75 converted to AB system by a 2.29 mag offset (Cowie, 
Songaila, & Barger 1999). As was seen in the rest B band LF, the rest 2000 A LFs also 
show a luminosity evolution that is stronger than in the rest B LF and continues up to 
2.5 < z < 3.5 bin. The M| 000 value is getting luminous from z = 0.6 to z — 3.5 over 1 mag, 
which is more significant than M| 000 errors. 

Cowie, Songaila, & Barger (1999) have derived the rest 2000A LF at 0.5 < z < 0.9 and 
1.0 < z < 1.5. We represent their LFs in the case of a = —1.0 and h = 0.75 as dashed line 
in corresponding redshift bins of Figure 7. They mentioned that their value of ultraviolet 
luminosity density, which is lower than that of other studies, would partly be caused by field 
variation. Taking account of this normalization fluctuation, as well as their fairly large error 
of ±0.4mag in determination of M*, our results agree with their LFs. 

For 2.5 < z < 3.5 in Figure 7, the dashed line represents the rest 1700A LF of Lyman- 
break galaxies at (z) = 3.04 by Steidel et al. (1999). We ignore the rest-band difference 
between 2000 and 1700A because the correction would be less than SM ~ 0.13. Again as 
seen in the rest B LF, our estimate of the rest UV LF at z ~ 3 is in good agreement with 
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that of Steidel et al. (1999), who have estimated only the bright end (Mi 70 o < —20). 

As has been noticed in rest B band LF, we can see less faint galaxies at higher redshift 
in this rest UV LF. This is probably caused by the fact that our JT'-limited sample is 
insensitive to blue dwarf galaxies. The deficit of faint-end galaxies in the if-selected sample 
is also seen in local universe. The LFs of if-selected sample (e.g., Cole et al. 2001; Table 3) 
have shallower faint-end slopes a ~ —1.0 than those of optical-selected sample (e.g., Blanton 
et al. 2001), a 1.20. 

5. Evolution of Rest XJV-K' Color 

In the previous section, concerned with rest-frame luminosity functions, we found the 
rest if'-band LF has almost no evolution, while the B- and the UV-band LFs show drastic 
luminosity evolution; that is, galaxies get brighter in bluer bands with increasing redshift. As 
a natural consequence, galaxies are inferred to be getting bluer in rest frame with redshift. 
To see this trend more definitely, we show the rest M2000 ~ Mk> color evolution with redshift 
in Figure 8. The rest-frame 2000A and K' luminosity of each galaxy selected in K' band can 
be calculated with the same manner described above. Note that this color estimate is not 
based on the same aperture magnitudes in 2000A and K' bands. We categorized the sample 
galaxies according to their absolute rest i^'-band magnitudes, M K >, which are fairly good 
indicators of the stellar mass. We have mentioned in §4.1 that M K > values derived from the 
extrapolation of SED template might be a poor proxy of stellar mass at z ~ 3, as suggested 
by Shapley et al. (2001) and Papovich, Dickinson, & Ferguson (2002). Though predominant 
uncertainty caused by the second old population would increase the absolute stellar mass 
estimate systematically, the correlation between Mk> and stellar mass would be maintained 
if all the galaxies at z ~ 3 are assumed to have similar mass ratios of the young to old stellar 
populations. It was found that the derived stellar masses have a fairly tight correlation 
with extinction-uncorrected rest V luminosity (Shapley et al. 2001). Papovich, Dickinson, 
& Ferguson (2002) have also found the correlation even with the rest UV luminosity with a 
scatter of ~ 0.3dex, though Shapley et al. (2001) have found only weak correlation. However, 
we should note again that the extinction is still uncertain in the UV luminosity at z ~ 3. 
In Figure 8, the most massive Mk> < —23.5 sample is represented by open circles. The 
—23.5 < Mx> < —21.0 and —21.0 < Mr-/ samples are represented by filled circles and 
crosses, respectively. The median M 2 ooo — Mk - ' color for each redshift is traced by solid, 
dotted and dashed lines from the most massive sample to the least massive one, respectively. 

Two significant trends can be seen in this figure. First, all three subsamples have a 
mild bluing trend in M 2 ooo ~ M K i with increasing redshift as expected from our results on 
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rest LFs. For example, M K > < —23.5 sample has median color M 2 ooo — M K / = 7.46 around 
z = 0.3, while at z — 2.7 there are no such red galaxies and median color becomes as blue 
as 4.73. The bluing trend we found has already been indicated by several studies (e.g., 
Cowie et al. 1996; Papovich, Dickinson, & Ferguson 2002) and is as qualitatively expected 
from a passive evolution model. Papovich, Dickinson, & Ferguson (2002) suggests that 
the UV-optical rest-frame SEDs of the Lyman break galaxies at 2.0 < z < 3.5 are much 
bluer than those of present-day spiral and elliptical galaxies. The other two subsamples, 
—23.5 < M K i < —21.0 and —21.0 < M K > also show this bluing trend, at least at z < 1.5, 
and almost no color changes have been found beyond z = 1.5. Our sample is so small that 
we cannot discriminate the difference of rapidity of this bluing trend among subsamples. 

It should be noted that an observational limit might affect this diagram. It is difficult to 
show explicitly the observational limit of K' < 24.0 in this rest-frame plane. For simplicity, 
we show in Figure 8 a fiducial line (heavy solid line) denoting the bluer color limit of a galaxy 
having constant luminosity of M2000 — — 20, which can be observed in the i^'-band image 
[K' < 24.0), ignoring the small k— correction (e.g., Mannucci et al. 2001) at the K' band. 
Galaxies brighter than M 2 ooo — ~ 20 would be detectable above this line. For a galaxy with 
^2000 = — 19, you may shift the lower limit upward by lmag. At higher redshift, it is found 
that galaxies that are bluer and fainter in 2000A are not likely to stay in this X'-limited 
sample. Our subsample of —23.5 < M K < —21.0 and —21.0 < M K would have been affected 
by this observational bias which may cause a lack of blue galaxies at high redshift if such 
blue population really exists. The blueing trend of these subsample seems to stop at z — 1.5 
in Figure 8, though it may be possibly attributed to this bias. 

On the other hand, for M K < —23.5 subsample, which seems to be less affected by this 
bias, all the galaxies become redder with time in z < 1, while there are galaxies around 
z ~ 1 as blue as those at z ~ 3. It seems that the scatter of the color becomes larger 
in the nearby universe. The same trend can be seen in the rest-frame U — B color of 
the brightest {Mb < —20.3) galaxy sample in HDF by Dickinson (2000). These galaxies 
with similar stellar masses keep their blue color throughout the passive evolution during a 
longer timescale than the typical timescale of an intense starburst event. This implies that a 
merging process (e.g., Lacey & Cole 1993) and subsequent star formation have been taking 
place continuously during 1 < z < 3 epoch. These blue galaxies of M K < —23.5 at z — 1 
seem to be assembled from lower mass galaxies via. merging and are ongoing star formation. 

Second, it is found that more massive galaxies tend to be redder in rest M 2 ooo — Mr 7 
color at every epoch. The median color of massive galaxies is always redder than that of less 
massive galaxies. The trend that red galaxies are more luminous than blue ones has already 
been established in the local universe for example by Blanton et al. (2001), who found that 
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on average high surface brightness, red, and highly concentrated galaxies are more luminous 
than low surface brightness, blue, and less concentrated galaxies. If rest M 2 ooo — M K i color 
indicates the age of a galaxy, the massive galaxies tend to have older, which can be easily 
inferred from the context of downsizing hypothesis, in which massive (spheroidal) galaxies 
formed at high z, followed at later epochs by the formation of less massive galaxies. 

6. Contribution to The Global Star Formation Rate Density 

In the last few years, there have been great advances in understanding the evolution of 
the global star formation rate density (SFRD). At low redshift between z = and z = 1.5, 
the SFRD is found to be steadily increasing with redshift. At high redshift, some results 
(e.g., Madau, Pozzetti, & Dickinson 1998) meet the SFRD peak around z ~ 1.5, while others 
(e.g., Steidel et al. 1999) maintain it relatively constant. The estimate of SFRD based on 
the ultraviolet luminosity of galaxies is strongly affected not only by dust obscuration, which 
many investigators took into account but also by the faint end of the LF. At last in this 
study, we derive the SFRD based on our evaluated rest UV LF. 

The SFRD can be estimated by integrating over the best-fit Schechter function to infinity 
for the rest UV LF in each redshift bin. Equation (4) below was used to convert luminosity 
at 2000A, 

L 2 ooo to an SFR based on the calibration given by Madau, Pozzetti, & Dickinson 

(1998) for a Salpeter IMF, 

SFRiMoyr- 1 ) = 1.26 x Kr 28 L 200 o(ergs s^Rz' 1 ). (4) 

The result is represented in Figure 9, in which several previous estimates are also plotted. 
All the previous results except Steidel et al. (1999) were derived by integrating the Schechter 
function to infinity. For consistency, we increased the Steidel et al. (1999) points by an 
additional factor 1.7 as they have noted in the case of integration to infinity. We checked 
our data point by recalculating the SFRD from other rest UV wavelengths such as 1500 and 
2800A as was adopted in previous studies and found that our SFRD estimate is stable within 
a fluctuation of 0.2. The dust extinction is uncorrected for in all the estimates. 

Our results are denoted as solid circles and are broadly consistent with others at z < 
1.5, though they show slight deficiencies at z > 1.5, where there are yet no established 
observations. This has probably arisen because our estimated rest UV LF of the ^'-selected 
sample has a shallower faint-end slope than that of Steidel et al. (1999). For an example, 
if we change the LF slope at 2.5 < z < 3.5 from a = —0.83 to a = —1.6 as Steidel et al. 

(1999) have adopted at z — 3, our estimate at z = 3 in Figure 9 shifts up to the open circle, 
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which is almost consistent with Steidel et al. (1999). Otherwise, if the slope of the rest UV 
LF at z = 3 is really as steep as a = —1.6, our K' < 24.0 limited sample contributes only 
42% of the total UV luminosity density, that is, of the total star formation rate density at 
z = 3. The remaining 58% should be charged to faint galaxies with M 2 ooo > — 19 for which 
our sample seems to have a deficiency of galaxies compared with the LF derived by Steidel et 
al. (1999). These missed faint galaxies correspond to very blue galaxies, which are probably 
affected by the observational limit in K' band as has been shown in §5. However, we should 
note here that as they have mentioned the a derived by Steidel et al. (1999) has a large 
uncertainty because it is determined only at the bright end (Mi 70 o < —20) of the LF. Poli 
et al. (2001) have estimated the photometric redshifts for /-selected galaxies and obtained a 
shallower faint end slope a = —1.37 ± 0.20 for the rest 1700A LF of 2.5 < z < 3.5 galaxies 
down to M 170 o = —19. The existence of a substantial population of such faint (M 170 o > —20) 
and blue Lyman break galaxies at z > 3 is still questionable. 

7. Summary 

We have constructed a very deep if'-selected multicolor BVRIz'JK' sample of field 
galaxies on the Subaru Deep Field. Based on this multicolor data, a photometric redshift 
for each sample galaxy was estimated. Our principal conclusions are the following: 

1. The overall N(z) distribution for the K' < 21 sample is consistent with previous obser- 
vations, and we for the first time derived the N(z) down to K' = 24. Taking account 
of dust extinction and selection effects of the sample, the observed N(z) distribution at 
the most secure magnitude bin, 22 < K' < 23, is well described with the PLE model, 
while a recent version of the hierarchical galaxy formation model shows an apparent 
deficiency of galaxies, especially at z > 2. The PLE model with realitic dust evaluation 
can account for the abundance of the high z (z > 2) galaxies. In the faintest bin of 
N(z), 23 < K' < 24, both of the model predictions are found to be inconsistent with 
the observations. The discrepancy in the faintest bin is presumably caused by the 
photometric errors that give an incorrect solution of photometric redshifts. 

2. The rest K', B, and UV(2000A) band luminosity functions and their evolution have 
been investigated. The rest K'-band luminosity function shows almost no evolution 
up to z = 3. The rest B luminosity function shows mild brightning with epoch up to 
z = 2.5, though its brightening degree is within the M* error caused by photometric 
and redshift errors. The rest 2000A LFs also show a luminosity evolution that is more 
significant than in the rest B LF and continues up to z = 3.5. The M| 000 value is 
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getting luminous from z = 0.6 to z = 3.5 over 1 mag, which is more significant than 
M 2*ooo errors - 

3. The evolutionary trend of rest luminosity functions appears to correspond to the bluing 
trend of the rest UV — K' color with increasing redshift. We also found that more 
massive galaxies are redder at every epoch. 

4. The findings about the evolution of rest luminosity functions are qualitatively in favor 
of the PLE model. The star formation activity in a galaxy changes the rest blue 
luminosity, while the amount of old stellar population does not change. The resulting 
rest .fT'-band luminosity function showing almost no evolution up to z = 2.5 can also 
be consistent with the hierarchical galaxy formation model if M/L K decreases with 
redshift; that is, the number density of K'-band luminosity-selected galaxies does not 
significantly change with redshift, while the number density of stellar mass-selected 
galaxies decreases with look-back time. 

5. We have evaluated the contribution to the cosmic star formation rate density of our 
if'-selected galaxies. The star formation rate density derived from integration of our 
rest UV luminosity function is about 42% of that derived from Lyman break galaxies 
at z — 3. The result suggests that a large fraction of global star formation rate density 
at z > 2 could represent the contribution of faint blue galaxy population, which has 
not yet explicitly been identified in a high z universe. 

Though the present results rely heavily on an availability of photometric redshifts, we 
are confident that we have successfully derived a coherent picture of evolution of K' , B, 
and UV rest-frame luminosity functions over the range 0.6 < z < 3.5 to the faint end 
below the limit of the today's spectroscopic studies. Our result raises the question of how 
such massive galaxies have formed at high redshift with efficient star formation, which is 
apparently consistent with the PLE picture in the K band, within the framework of the 
structure formation in the CDM universe. Of course, deeper NIR imaging surveys in future 
will allow us to derive more secure results for N(z) at the faintest magnitude and on the faint 
end of the rest LFs, in which we have obtained robust results in this study. Deep observation 
at longer wavelength with SIRTF and ASTRO-F will be also required to constraint the 
SEDs of galaxies at z > 3. The technique that we have adopted in this study will be applied 
to future wide-field i^-band surveys to determine the rest luminosity at any wavelength at 
any redshift down to the faintest magnitude. On the other hand, ongoing large multiobject 
spectroscopic surveys are of course extremely important to reveal the galaxy evolution more 
reliably. 
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Fig. 1. — Comparison between z spec and z p hot for our spectroscopic sample of SDF. 
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Fig. 2. — K 1 magnitude vs. z phot diagram. The error bars of photometric redshifts are 
evaluated from the x 2 probability distribution at confidence level 90%. 




Fig. 3. — Redshift distribution of our SDF sample with different limiting magnitudes of 
K' < 23.0, K' < 23.5, and K' < 24.0. 
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Fig. 4. — Top to bottom: Normalized cumulative redshift distribution for our SDF sample 
for the bright-magnitude bin to the faint magnitude-bin. The histograms show the results of 
observation. The shaded regions show the ±3<r deviated counts estimated by Monte Carlo 
realizations when photometric redshift errors are taken into account. The heavy dashed 
lines denote the prediction of the PLE model by Totani et al. (2001b), while the heavy 
dotted lines denote the the hierarchical galaxy formation model by Nagashima et al.(2003). 
Top: Prediction of PLE and hierarchical clustering models in KC98 as well. 
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Fig. 5. — Top to bottom: Rest if'-band luminosity functions at z — 0.6 — 1.0, z = 1.0 — 1.5, 
z = 1.5 — 2.5, and z = 2.5 — 3.5. Filled circles show the luminosity function determined in 
the present study, and solid lines are the best-fit Schechter functions. The dotted line shows 
the local K' luminosity function by Loveday (2000) 
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Fig. 6. — Same as Figure 5 for rest I?-band LFs. The local i?-band luminosity function 
denoted as dotted line is from Blanton et al. (2001), converted from their g* band to B 
band. Top: 0.75 < z < 1.0 rest 5-band LF derived on CFRS (Lilly et al. 1995) is also 
superimposed as a dashed line. Bottom: In the 2.5 < z < 3.5 bin, a dashed line shows the 
rest V^-band LF for Lyman-break galaxies at (z) = 3.04 by Shapley et al. (2001), converted 
to rest-B band and fl m = 1.0, h = 0.75 cosmology. 
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Fig. 7.— Same as Figure 5 for rest UV(2000A AB)-band LFs. The local 2000A luminosity 
function by Sullivan et al. (2000) for h = 0.75 is denoted as dotted line converted from their 
magnitude system to AB system with a 2.29 mag offset. In the bins for 0.6 < z < 1.0 and 
1.0 < z < 1.5, dashed lines show the rest 2000A LF for 0.5 < z < 0.9 and 1.0 < z < 1.5 
derived by Cowie, Songaila, & Barger (1999) for a = —1.0 and h = 0.75. In the bin for 
2.5 < z < 3.5 panel, the dashed line shows the rest 1700A LF for Lyman-break galaxies at 
(z) = 3.04 by Steidel et al. (1999). 
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Fig. 8. — Rest M 20 oo — Mk> color evolution with redshift. Open circles, filled circles and 
crosses denote the subsamples with respect to the absolute K' magnitudes, M K i < —23.5, 
—23.5 < Mk> < —21.0 and —21.0 < Mk>, respectively. The median M2000 — Mk> color 
for each redshift is traced by solid, dotted and dashed lines from bright subsample to faint 
subsample, respectively. We show a thick solid line as a fiducial line denoting the lower 
color limit corresponding to our sample's magnitude limit of K' lim = 24.0 galaxy for constant 
-^2000 = — 20 ignoring the small k— correction at f^'-band. Galaxies brighter than M 20 oo 

20 would be detectable above this line. If M2000 = — 19 galaxy, you may shift the lower 
limit upper by lmag. 
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Fig. 9. — Star formation rate density as a function of redshift (Madau plot) derived from the 
integration of UV-LF. The estimates of this work are shown in filled circles. The different 
data are taken from previous estimates of Lilly et al. (1996; open squares), Connoly et al. 
(1997; crosses), Cowie et al. (1999; filled triangles), Madau et al. (1998; open triangles), 
and Steidel et al. (1999; filled squares). The dust extinction is uncorrected for in all the 
estimates. If we change our rest UV LF slope from a = —0.83 to a = —1.6 as Steidel et al. 
(1999) have adopted at z = 3, the estimate at z = 3 shifts up to the open circle from the 
filled circle. 
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Table 1. Summary of Optical Imaging Data for SDF. 



band 


exp.time 


seeing size 






(sec) 


(arcsec) 


(mag) 


B 


9000 


0.6 


28.08 


V 


3600 


0.8 


26.28 


R 


14400 


0.7 


27.88 


I 


3600 


0.7 


25.93 


z' 


3600 


0.8 


24.27 



a The 5<T limiting magnitudes in the Vega 
system within 2" aperture 



Table 2. Chance Probability and \ 2 ( m parentheses) Values Obtained from 
Kolmogorov-Smirnov Test for N(z) Distribution. 



models 


19 < K' < 21 


21 < K' < 22 


22 < K' < 23 


23 < K' < 24 


PLE 


1.97%(8.90) 


1.03%(10.2) 


18.2%(4.67) 


0.232%(13.3) 


hierarchical 


9.32%(5.91) 


76.4%(1.72) 


6.46e-4%(24.7) 


4.49e-13%(66.2) 
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Table 3. The Best-Fit Shechter Parameters of Rest-Frame Luminosity Functions in Each 

Redshift Bin Derived on SDF 



Rest-frame z range a M* <p* 



K' 


0.6- 
1.0- 


1.0 
1.5 


i qcr+0.05 
i - OJ +0.02 
_1 oc-0.01 
i -°°-0.10 


-25.03+°^ 
-25.02185 


O 07+2.12 

d -°'-0.41 

q O7 + 0.94 

d -°'-2.18 


x 10" 
x 10" 


-3 
-3 




1.5- 


2.5 


07+O.O9 
J-- ' +0.00 


-24.95ftS 


q 1 q+2.09 

°- io +0.09 


x 10- 


-3 




2.5 - 


3.5 


-1.70_ 006 


-23.951°;- 


2-5711;?! 


x 10- 


-3 


B 


0.6- 


1.0 


-1 07+ - 01 

L - u 1 -0.07 


-20.2618;^ 


i-76:8:^ 


x 10" 


-2 




1.0- 


1.5 


-0 Q2+ - 13 
u - yz -0.10 


-2o.6it8: 3 ? 


1.8018; 3 7 8 


x 10- 


-2 




1.5- 
2.5- 


2.5 
3.5 


_1 1O+0.25 

1 - ±a +o.oi 
-1.49t°:£ 


-2o.82±8:SS 


O qi +7.82 
°- 0± +1.76 

3.86±l:2S 


x 10" 
x 10' 


-3 
-3 


[/(2000A) 


0.6- 


1.0 


-0.9918;?* 


-18.49l8;n 


2 70+0.80 

z - ' a -0.08 


x 10" 


-2 




1.0- 


1.5 


-o.75±8:S2 


i o 7i +0^07 


3.49:8;?* 


x 10- 


-2 




1.5- 


2.5 


-0-83±8;S 


-19.1218;^ 


1 qc;+o.88 


x 10- 


-2 




2.5- 


3.5 


-o.83ig:S 


-19.5418;^ 


i.nl8;8I x 10- 


-2 



